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Abstract 

Background: Blood-sucking lice (suborder Anoplura) parasitize eutherian mammals with 67% of the 540 described 
species found on rodents. The five species of blood-sucking lice that infest humans and pigs have fragmented 
mitochondrial genomes and differ substantially in the extent of fragmentation. To understand whether, or not, 
any life-history factors are linked to such variation, we sequenced the mt genomes of Polyplax asiatica and Polyplax 
spinulosa, collected from the greater bandicoot rat, Bandicota indica, and the Asian house rat, Rattus tanezumi, 
respectively. 

Results: We identified all of the 37 mitochondrial genes common to animals in Polyplax asiatica and Polyplax 
spinulosa. The mitochondrial genes of these two rat lice are on 1 1 circular minichromosomes; each 
minichromosome is 2-4 kb long and has 2-7 genes. The two rat lice share the same pattern for the distribution of 
the protein-coding genes and ribosomal RNA genes over the minichromosomes, but differ in the pattern for the 
distribution of 8 of the 22 transfer RNA genes. The mitochondrial genomes of the Polyplax rat lice have 3.4 genes, 
on average, on each minichromosome and, thus, are less fragmented than those of the human lice (2.1 and 2.4 
genes per minichromosome), but are more fragmented than those of the pig lice (4.1 genes per minichromosome). 

Conclusions: Our results revealed distinct patterns of mitochondrial genome fragmentation within the genus 
Polyplax and, furthermore, indicated a possible inverse link between the extent of mitochondrial genome 
fragmentation and the length of life cycle of the blood-sucking lice. 
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Background 

Lice in the suborder Anoplura are wingless, exclusive 
blood-sucking insects and are permanent ectoparasites 
of eutherian mammals [1,2]. Blood-sucking lice evolved 
from chewing lice -100 million years ago (Mya) and 
diversified rapidly -65 Mya with their mammalian hosts 
[3]. More than 540 species of blood-sucking lice have 
been described and are classified into 15 families 



* Correspondence: daochaojin@126.com 

institute of Entomology, Guizhou University, and the Provincial Key 
Laboratory for Agricultural Pest Management in Mountainous Region, 
Guiyang 550025, China 

Full list of author information is available at the end of the article 

(3 BioMed Central 



and 50 genera [4-7]. Twelve of the 29 recognized 
mammalian orders and -840 mammalian species are 
hosts of blood-sucking lice. The diversity of blood-sucking 
lice mirrors that of their mammalian hosts with 67% of 
the described species found on rodents [2,5,7]. Each 
species of mammalian host is usually parasitized by a 
single species of blood-sucking lice, but there are 
exceptions as many mammalian species are hosts of 
multiple species of blood-sucking lice (up to seven 
species) [7]. Vice versa, each species of blood-sucking 
lice usually parasitizes only one species of mammals, but 
many species of blood-sucking lice parasitize multiple 
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species of mammals that are closely related phylogenetically 
(up to 39 species) [6]. 

Adding to their unique, exclusive blood-feeding life style 
and the high host-specificity (relative to other ectoparasites), 
recent studies revealed that blood-sucking lice have an un- 
usual, fragmented mitochondrial (mt) genome organization. 
For insects and other bilateral animals, the mt genome is 
usually a single, circular chromosome, 13-20 kb long, with 
37 genes: 13 for proteins, two for ribosomal RNAs and 22 
for transfer RNAs [8,9]. The mt genomes of the human 
body louse, Pediculus humanus, and the human head louse, 
Pediculus capitis, however, have 20 minichromosomes; each 
minichromosome has 1-3 genes and is 3-4 kb in size 
[10,11]. The 34 mt genes identified in the human pubic 
louse, Pthirus pubis, are on 14 minichromosomes; each 
minichromosome has 1-5 genes and is 1.8-2.7 kb in size 
[10,11]. The domestic pig louse, Haematopinus suis, and 
the wild pig louse, Haematopinus apri, have their 37 mt 
genes on nine minichromosomes; each minichromosome 
has 2-8 genes and is 3-4 kb in size [12]. 

The fragmented mitochondrial genomes of the human lice 
and the pig lice represent the most radical departure to date 
in bilateral animals from the typical, single-chromosome 
organization of mt genomes, although multipartite mt 
genomes have also been observed in the rotifer, Brachionus 
plicatilis, and the booklouse, Liposcelis bostrychophila, which 
have two mt chromosomes [13,14], and the potato cyst 
nematode, Globodera pallida, which has six chromosomes 
[15]. Outside bilateral animals, highly fragmented mt 
genomes were observed in ichthyosporean protists 
[16], diplonemid protists [17] and box jellyfish [18]. 

Although all having fragmented mt genomes, the 
human lice and the pig lice differ substantially in the 
extent of mt genome fragmentation and in the distribution 
of mt genes over the minichromosomes [12]. To understand 
whether, or not, any life-history factors are linked to such 
variation, we sequenced the mt genomes of two species of 
blood-sucking lice that parasitize rats, Polyplax asiatica 
(Ferris 1923) and Polyplax spinulosa (Burmeister 1839). 
Polyplax is one of the most species-rich genus of the 
suborder Anoplura with 78 described species, second only 
to the genus Hoplopleura, which has 141 described species 
[6]. Polyplax species shared their most recent common 
ancestor (MRCA) with primate lice -47 Mya and their 
MRCA with pig lice -67 Mya [3]. Po. asiatica infests the 
greater bandicoot rat, Bandicota indica, three other 
related species of rats and occasionally the Asian house 
shrew, Suncus murinus [7,19]. Po. spinulosa infests the 
Asian house rat, Rattus tanezumi, the brown rat, Rattus 
norvegicus, the black rat, Rattus rattus, and six other 
related species of rats [6,7]. While Po. asiatica is only 
found in Asia, Po. spinulosa has a worldwide distribution 
and is responsible, as a vector, for the transmission of patho- 
genic microorganisms such as Mycoplasma haemomuris 



(formerly known as Haemobartonella muris), Rickettsia 
typhi, Trypanosoma lewisi, Borrellia duttoni and Brucella 
brucei [20-22]. We show that both Po. asiatica and Po. 
spinulosa have fragmented mt genomes; these two 
species, however, differ in the pattern of mt genome 
fragmentation. Furthermore, our comparison among the 
rat lice, the human lice and the pig lice indicates a 
possible inverse link between the extent of mt genome 
fragmentation and the length of life cycle of these lice. 

Methods 

Collection of rats and lice 

The blood-sucking lice, Po. asiatica and Po. spinulosa, 
were collected in 2011 from the greater bandicoot 
rat, B. indica, and the Asian house rat, R. tanezumi, 
respectively, in Jinping county, Yunnan province, China. 
The rats were caught with trap-cages set indoors (farmers' 
houses, barns and stables) and outdoors (farmlands, 
scrublands and woodlands). Alive rats trapped were placed 
individually in pre-marked cotton bags and transferred to 
laboratory for species identification and parasitological 
check. Blood-sucking lice on the body surface of each 
rodent host were collected and preserved in 95% ethanol 
at -20°C prior to DNA extraction. Samples of Po. asiatica 
and Po. spinulosa and their rat hosts were deposited in the 
Institute of Pathogens and Vectors, Dali University. The 
capture of rodents was approved by health authorities in 
Yunnan province, China. Animal protocols and procedures 
were approved by the animal ethics committees at 
Guizhou University and Dali University. 

DNA extraction, mitochondrial genome amplification and 
sequencing 

Genomic DNA was extracted from individual louse 
specimens with DNeasy Tissue kit (QIAGEN). For Po. 
asiatica (sample #57) collected from the greater bandicoot 
rats, a 532-bp fragment of coxl and a 360-bp fragment of 
rrnL were initially amplified by PCR with primer pairs 
mtd6-mtd9 andl6SF-Lxl6SR (see Additional file 1). These 
two pairs of primers target conserved sequence motifs in 
coxl and rrnL; the amplicons were sequenced directly 
using Sanger method at the Tiangen Biotech, Beijing 
(TBB). Two pairs of specific primers for Po. asiatica, 
cox57F-cox57R and 16S57F-16S57R, were designed from 
sequences of the coxl and rrnL fragments. The two 
specific primers in each pair go outbound with 
131 bp and 1 bp, respectively, in between. PCRs with 
these specific primers amplified two near full-length 
mt minichromosomes of Po. asiatica that contain 
coxl and rrnL respectively; these amplicons (3.2 kb and 
2.6 kb in size) were sequenced using Sanger method at 
the TBB. Another pair of primers specific to Po. asiatica, 
57F-57R, was designed from conserved non-coding 
sequences that flank the coding regions of the two 
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minichromosomes above. The PCR with 57F-57R primer 
pair produced a mixture of amplicons ranging from 1 to 
2 kb in size, expected from the coding regions of all mt 
minichromosomes of Po. asiatica (Figure 1A). These 
amplicons were sequenced with Illumina Hiseq 2000 
platform at the Beijing Genome Institute, Hong Kong 
(BGI-HK). The PCR strategy used in this study was 
developed from the observations we made in previous 
studies on the human lice and the pig lice that each 
mt minichromosome has a distinct coding region but 
a well-conserved non-coding region [10-12]. 

For Po. spinulosa (sample #301) collected from the 
Asian house rats, a 452-bp fragment of rrnS and a 
360-bp fragment of rrnL were amplified initially by 
PCR with primer pairs 12SA-12SB and 16SF-Lxl6SR 
(see Additional file 1). These two pairs of primers target 
conserved sequence motifs in rrnS and rrnL and the PCR 
amplicons were sequenced directly using Sanger method 
at the TBB. Two pairs of specific primers for Po. spinulosa, 
12S301F-12S301R and 16S301F-16S301R, were designed 
from sequences of the rrnS and rrnL fragments. The two 
specific primers in each pair go outbound and are 1 bp 
and 89 bp respectively from each other. PCRs with 
these two pairs of specific primers amplified two mt 
minichromosomes of Po. spinulosa that contain rrnS and 
rrnL respectively; these amplicons (2 kb and 2.3 kb in size) 
were sequenced using Sanger method at the TBB. Another 
pair of primers specific to Po. spinulosa, 301F-301R, was 
designed from conserved non-coding sequences that flank 
the coding regions of the two minichromosomes above. 
The PCR with primer pair 301F-301R produced a mixture 
of amplicons ranging from 1 to 2 kb in size, expected from 
the coding regions of all mt minichromosomes of Po. 
spinulosa (Figure 1A). These amplicons were sequenced 
with Illumina Hiseq 2000 platform at the BGI-HK. 



Taq DNA Polymerase (Tiangen Biotech) was used in the 
initial short PCRs with the following cycling conditions: 
94°C for 1 min; 35 cycles of 98°C for 10 sec, 45°C for 
30 sec, 72°C for 1 min; and a final extension of 72°C for 
2 min. LA Taq (TaKaRa) was used in the long PCRs with 
the cycling conditions: 94°C for 1 min, 35 cycles of 98°C 
for 10 sec, 60-65°C (depending on primers) for 30-40 sec, 
68°C for 3 min; and 72°C for 6 min. Positive and negative 
controls were executed with each PCR. PCR amplicons 
were checked by agarose gel (1%) electrophoresis; the sizes 
of PCR amplicons were estimated by comparing with 
molecular markers. PCR products were purified with 
Wizard SV Gel and PCR clean-up system (Promega). 

Assembly of Illumina sequence-reads, gene identification 
and verification of individual mitochondrial 
minichromosomes 

Illumina sequence-reads obtained from the mt minichro- 
mosomes of Po. asiatica and Po. spinulosa were assembled 
into contigs with Geneious 6.1.6 [23]; the parameters for 
assembly were minimum overlap identity 98% and 
minimum overlap 50 bp. tRNA genes were identified 
using tRNAscan-SE [24] and ARWEN [25]. Protein-coding 
genes and rRNA genes were identified with BLAST 
searches of GenBank [26,27]. Identical sequences shared 
between genes were identified with Wordmatch [28]. 
Sequence alignments were with Clustal X [29]. The size 
and circular organization of each mt minichromosome of 
Po. asiatica and Po. spinulosa identified by sequence-read 
assembly were verified by PCR (Figure IB, 1C). A pair of 
outbound primers (forward and reverse) was designed 
from the coding region of each minichromosome 
(see Additional file 2). The two primers in each pair 
were next to each other with a small gap or no gap in 
between. PCRs with these primers amplify the full or near 
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Figure 1 PCR amplification of the mitochondrial (mt) minichromosomes of the Polyplax rat lice. (A) Lane 1: GeneRuler®100 bp DNA 
Ladder (Thermo Scientific). Lane 2: PCR amplicons generated with primer pair 57F-57R that spans the coding region of each mitochondrial (mt) 
minichromosome of Polyplax asiatica. Lane 4: PCR amplicons generated with primer pair 301 F-301 R that spans the coding region of each mt 
minichromosome of Polyplax spinulosa. (B) PCR verification of the mt minichromosomes of Po. asiatica. Lane 1 and 13: 500 bp DNA Ladder 
(JIANGEN). Lane 2-12: PCR amplicons from the 1 1 minichromosomes of Po. asiatica: atp8-atp6, trnE-cob-trnl, cox1-tmL 2 (taa), trnD-tmY-cox2-nad6, 
trnR-nad4L-cox3-trnA, trnS 1 (tct)-trnS 2 (tga)-nadl-tmT-tmG-nad3-tmW, trnQ-nad2-tmN-tmP, tmK-nad4-trnF, trnH-nad5, rrnS-trnC, trnM-trnL^tagyrrnL-trnV. 
(C) PCR verification of the mt minichromosomes of Po. spinulosa. Lane 1 and 13: 500 bp DNA Ladder (JIANGEN). Lane 2-12: PCR amplicons from the 1 1 
minichromosomes of Po. spinulosa: atp8-atp6, trnE-cob-trnl, coxl-trnL^tag), trnT-trnD-trnY-cox2-nad6-trnA l trnR-nad4L-trnP-cox3, nadl-trnG-nad3-trnW, 
trnQ-nad2-trnN, trnK-nad4, trnH-nad5-trnF, trnS 1 {tct)-trnS 2 {tga)-rrnS-trnC, trnM-trnL2{taa)-rrnL-trnV. Genes from which PCR primers were designed 
are in bold. 



Dong et al. BMC Genomics 2014, 15:44 
http://www.biomedcentral.eom/1 471 -21 64/1 5/44 



Page 4 of 12 



full length of each minichromosome if it has a circular 
organization. PCR set-up, cycling conditions, agarose gel 
electrophoresis and size measurement were the same 
as described above. The nucleotide sequences of the mt 
genomes of Po. asiatica and Po. spinulosa have been de- 
posited in GenBank under accession numbers KF647751- 
KF647772. 

Results 

Mitochondrial genome of Polyplax asiatica, the louse of 
the greater bandicoot rats 

We obtained 546,066 sequence-reads from the mt 
genome of Po. asiatica by Illumina sequencing (Table 1). 
The sequence-reads are 90 bp each; assembly of these 
sequence-reads into contigs allowed us to identify all of 
the 37 mt genes typical of animals in Po. asiatica distributed 
over 11 minichromosomes (Figure 2 A; Table 1). Each 
minichromosome is circular, 2.4-3.3 kb in size (Figure IB), 
and consists of a coding region and a non-coding region. 
Each coding region contains two to seven genes, and 
varies in size from 802 bp for rrnS-trnC minichromosome 
to 1,756 bp for trnH-nadS minichromosome (Table 1; 
Figure 2A) (Note: minichromosomes are named after 
their genes hereafter). Seven of the 11 minichromosomes 
of Po. asiatica have one protein-coding or rRNA gene 
each; the other four minichromosomes have two protein- 
coding genes each. The 22 tRNA genes are scattered over 
10 of the 11 minichromosomes; each minichromosome 
has one to five tRNA genes except atp8-atp6, which 
has no tRNA genes (Figure 2A). Each of the 37 mt 
genes identified in Po. asiatica is present on only one 
minichromosome; there is no overlap in gene content 

Table 1 Mitochondrial minichromosomes of Polyplax 
asiatica and Polyplax spinulosa identified by Illumina 
sequencing 

Minichromosome Size of coding Number of Illumina 





region (bp) 


sequence-reads 


otp8-otp6 {otp8-otp6) 


838 (832) 


73276 (143560) 


E-cob-l (E-cob-l) 


1242 (1228) 


47008 (50108) 


cox1-L 2 {cox1-L 1 ) 


1626 (1599) 


49653 (49675) 


D-Y-cox2-nad6 {T-D-Y-cox2-nod6-A) 


1292 (1399) 


35911 (41691) 


R-nad4L-cox3-A {R-nad4L-P-cox3) 


1268 (1190) 


36763 (64203) 


S r S 2 -nad1-T-G-nad3-W 
{nad1-G-nad3-W) 


1602 (1419) 


21705 (66380) 


Q-nad2-N-P {Q-nad2-N) 


1385 (1134) 


24989 (46412) 


K-nad4-F {trnK-nod4) 


1465 (1331) 


38589 (78061) 


H-nad5 {trnH-nad5-trnF) 


1756 (1816) 


39554 (44417) 


rrnS-C (S, -S 2 -rrnS -Q 


802 (904) 


121142 (165740) 


M-L r rrnL-V (M-L 2 -rrnL-V) 


1333 (1237) 


57476 (96783) 


Total 


14609 (14089) 


546066 (847030) 



Note: Gene arrangements and numbers outside brackets are for Polyplax 
asiatica and those in brackets are for Polyplax spinulosa. 



between different minichromosomes. All of the 37 mt 
genes have the same orientation of transcription relative 
to the non-coding region, except trnT and nadl, which 
have the opposite orientation of transcription relative to 
other genes (Figure 2A). 

We sequenced the full-length non-coding regions, 
1,532 bp and 1,244 bp respectively, of two mt minichromo- 
somes of Po. asiatica-. coxl-trnL 2 and trnM-trnLj-rrnL-trnV. 
The length variation between the two non-coding regions is 
due to a 304-bp section that is present only in coxl-trnL 2 
minichromosome lying downstream the coding region and 
a 25-bp section that is present only in trnM-trnL r rrnL~ 
trnV minichromosome lying upstream the coding region 
(Figure 3A). Excluding these two sections, the non-coding 
regions of the two minichromosomes have 96% identity to 
each other. As in the human lice and the pig lice [11,12], an 
AT-rich motif (140 bp, 64% A and T) is present in the 
non-coding regions of Po. asiatica upstream the 5'-end of 
the coding region, whereas a GC-rich motif (76 bp, 71% G 
and C) is present downstream the 3'-end of the coding re- 
gion (Figure 3A). In addition to the full-length non-coding 
region sequences of the two minichromosomes, we also 
obtained partial sequences of the non-coding regions, 
101-169 bp and 49-449 bp respectively, upstream and 
downstream the coding regions of the other nine 
minichromosomes of Po. asiatica. Two highly conserved 
sequence-motifs, 101 bp and 44 bp long respectively, are 
present in the sections of the non-coding regions 
upstream and downstream the coding regions of all of 
the 11 minichromosomes (see Additional file 3). 

Mitochondrial genome of Polyplax spinulosa, the louse of 
the Asian house rats 

We obtained 847,030 sequence-reads from the mt genome 
of Po. spinulosa by Illumina sequencing (Table 1). As above 
for Po. asiatica, these sequence-reads are 90 bp each in 
length. We assembled these sequence-reads into contigs and 
identified all of the 37 mt genes typical of animals in Po. 
spinulosa. These genes are on 11 minichromosomes; each 
minichromosome is 2.0-2.9 kb in size and has a circular 
organization (Figure 1C; Figure 2B). As in Po. asiatica, seven 
of the 11 minichromosomes of Po. spinulosa have one 
protein- or rRNA-coding gene each; the other four mini- 
chromosomes have two protein-coding genes each. The 22 
tRNA genes are on 10 of the 11 minichromosomes; each 
minichromosome has one to four tRNA genes except 
atp8-atp6, which has no tRNA genes (Figure 2B). Each 
minichromosome has a coding region and a non-coding re- 
gion. The coding region of each minichromosome contains 
two to six genes, and varies in size from 832 bp for atp8- 
atp6 minichromosome to 1,816 bp for trnH-nad5-trnF 
minichromosome. With the only exception of nadl, all of 
the mt genes of Po. spinulosa have the same orientation of 
transcription relative to the non-coding regions (Figure 2B). 
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(A) Polyplax asiatica (louse of the greater bandicoot rat, Bandicota indica) 




(B) Polyplax spinulosa (louse of the Asian house rat, Rattus tanezumi) 




Figure 2 The mitochondrial (mt) genomes of Polyplax asiatica (A) and Polyplax spinulosa (B). Each minichromosome has a coding region 
(with gene name, transcription orientation and length indicated) and a non-coding region (in black). Minichromosomes are in alphabetical order 
by the names of their protein coding and rRNA genes. Abbreviations of gene names are: otp6 and otp8 (for ATP synthase subunits 6 and 8), 
coxl-3 (for cytochrome c oxidase subunits 1-3), cob (for cytochrome b), nadl-4 and nad4L (for NADH dehydrogenase subunits 1-6 and 4 L), rrnS 
and rrnL (for small and large subunits of ribosomal RNA). tRNA genes are shown with the single-letter abbreviations of their corresponding amino acids. 
Minichromosomes that have identical gene content and gene arrangement between the two Polyplax species are indicated with asterisk symbols "*". 



We sequenced the full-length non-coding regions, 
1,097 bp and 1,146 bp respectively, of two mt minichro- 
mosomes of Po. spinulosa: trnM-trnL 2 -rrnL-trnV and 
trnS r trnS 2 -rrnS-trnC. The length variation between the 
two non-coding regions is due to a 53-bp section that is 
present only in trnS r trnS 2 -rrnS-trnC minichromosome 



lying upstream the coding region (Figure 3B). Excluding 
the 53-bp section, the non-coding regions of the two 
minichromosomes have 98% identity to each other. 
As above in Po, asiatica, an AT- rich motif (192 bp, 
64% A and T) is present in the non-coding regions in 
Po. spinulosa upstream the 5'-end of the coding 
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0 
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GC-rich (65 bp, 57% C A C) 
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1 Adl IAAI ddl CI (AG! ICCICAdCCAI I AT IAI I CA( AACAI IAI AA IAAAA I AAA I AA ICAIAdl Al Cd(A(AI IAAAA 1 AA 1 AA 1 IAI AAA 1 A(A( Al AA 1 AAI AA 1 AAACA1 Al A( Al A( A(A(A( Adl Al A( ACAAM AACU A(. 


Al AA(.I A( AAA 1 Al (Al AA IAAAA lAlAIAAAA (Al Al AAA 1 dCAA 1 CI AACAAdC(AA(AA I Adl A 
AlAAdlACAAA IAI (AIAA IAAAA 1 A( A IAAAA (Al Al AAA 1 dCAA 1 C(AACAAdC( AA( AA 1 AdM 


A< AAAAA M(AAIAIC(.K.(A(.AAACCI IAI Al AA<A(AdCA(A1 (Al dl 1 1 Al A 1 ( Ad('l dCI Al A/ 
A< AAAAA (.l(AAIAIC(.l(.<A(AAACCI IAI AIAAIA(AdCA(AI (Al dl 1 IAI AlCAdCI <.CI Al A/ 


A IAIA( AC( A IddCd 
MIAIA(ACGATGGCG 


AlCMCAddddAA ICI (ACAICfAAdl (A(Adl CA( AAdd(A<Adl dl 1 (AA<ACACI CAdd<il AAA 
Al ('MCA(.d(i(.AA l(J < ACAI CdAA(.ICA(A(.l CA(AA<>d(A<AdI dl KAAIACACTCAGGGIAAA 


AA GACTCKAI ( A* Al Al Al Al 1 1 IAAACTCACI f GTACiTAACTCdl I AA 1 AIAA IAI AC(.( AAAAAA 
AA (iACICTGATGA(AIAIAI Al 1 1 IAAACTCACI T(iIAdTAACTC(iT TAAIATAAIAIACdCAAAAAA 


(AAA 1 C< AA 1 1 AAA 1 
(AAA 1 C< AA 1 1 AAA 1 


IAI AAI IAI AAA IAAI AAA 1 IAI AAAA lAdlAAdACAICICUAdCAAA IAICAI IAAAA 1 IAAI AA 
IAI AAI 1 Al AAA 1 AA 1 AAA 1 IAI AAAA 1 Adl AA<ACAI CI CMAdCAAA IAICAI IAAAA 1 1 AA 1 AA 


AA IAI ACAI AAAA (.(.(AC( AAAA 1 (.1 AAA 1 1 AACCC(A(A(,I lAdAA (.(AAAAA 1 KAAddll 1 IAI (.(AAA CCI IAAI AdAd 
AA 1 Al ACAI AAAA (.ddACIAAAA Idl AAA I IAACCC<A(AM 1 AIAA (.(AAAAA I KAAddl 1 1 lAIUAAACCI IAI IA(A(. 
AT-rich (192 bp, 64% A A T) 





AAIAI (Al AAAA ( Al CI I CI A I I IAI A(AC(AI d(<C(AICdl CAdddl (AAA I CI (ACAI C(AAdlCA( AM CI IAAAA (AAA dl dCAlAAddCACI 1 1 AdddCI I AAdl CCCI (Al (A(AI Al A(AI I ICAAdCICACI IdlddlAI I I 1047 
AA( Al ( Al AAAA CAI CI I t:i Al I IAI A( AC(A1 (.(.C( AlCdICAdddK AAA I CI (ACAI C(AAdTCA(AdICI IAAAA CAAA dl dCACAAddCACI I TAdddCI lAAdlCCCI (Al (A(AI Al A(AI I ICAA<>C! CACI I (.1 (.(.I A I I I 1045 

301 F ► 



Figure 3 Alignment of nucleotide sequences in the non-coding regions of the mitochondrial minichromosomes of Polyplax asiatica 
(A) and Polyplax spinulosa (B). 57F and 57R are the primers used to amplify the coding regions of all mitochondrial minichromosomes of Po. 
asiatica. 301 F and 301 R are the primers used to amplify the coding regions of all mitochondrial minichromosomes of Po. spinulosa. 



region, and a GC-rich motif (65 bp, 57% G and C) 
is present downstream the 3'-end of the coding region 
(Figure 3B). Additional to the full-length non-coding 
region sequences of the two minichromosomes, we 
also obtained partial sequences of the non-coding regions, 
95-217 bp and 74-136 bp respectively, upstream and 
downstream the coding regions of the other nine mini- 
chromosomes of Po. spinulosa. As in Po. asiatica, two 
highly conserved sequence-motifs, 95 bp and 74 bp long 
respectively, are present in the sections of the non-coding 
regions upstream and downstream the coding regions of 
all of the 11 minichromosomes (see Additional file 4). 



Discussion 

Variation in the pattern of mt genome fragmentation 
between the two Polyplax rat lice 

Only two of the 11 minichromosomes have identical gene 
content and gene arrangement between Po. asiatica and Po. 
spinulosa, revealing, for the first time, marked intra-genus 
variation in the pattern of mt genome fragmentation among 
blood-sucking lice (Figure 2). Comparison between the two 
Polyplax rat lice, and with the pig lice and the human lice, 
revealed the translocations of tRNA genes between mini- 
chromosomes in both Po. asiatica and Po. spinulosa after 
these two species split from their MRCA, and furthermore, 
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Species of lice 

Polypi® asiaka (louse 
Polyplax spinulosa (louse of Asian house ral 
Peiiculus knanus (human body louse) 
_| L Peiiculus capitis (human head louse) 
— Pthirus pubis (human pubic louse) 
r Mmtopiimapri (wild pig louse) 
L Haematopinus suis (domestic pig louse) 
Hypotheticalancestor of insects 



rC 



Gene content & gene arrangement in each mitochondrial minichromosome 

gSM-i, S/igHIW, fljpj-fl^, M-coxhai6, Mai4L^Lnai0M^^, 0-miWMC\ 
pxf : U E-cob-I, rf-Gwiarfj-II./ ^^ Mh-rrnkl 0-nc£-\,$i-5>rrnK. 



coil, cob, Q-nail, naiS, atpktpi lcoxl hxU\ l-mii K-nai{ lm-rml, F-naid, M, P-naiU, G-nai4L-Y,Li/2-rrn$-C, Si-M, T-M, ¥-k 
coxl, Mi, nail-Q, nail, atphtp6, J-cox2, \coxU\ G-naii-Y-K2, h-rrnl, F-nai6-E-M, P-nail-l, W-H-R-naiil, L-rrnS, C, 



\coxlM cox2I-D-atphtpkoxl-G-ndM4M 



Figure 4 Translocations of tRNA genes between mitochondrial minichromosomes and identity-swap between tRNA-leucine genes in 
the rat lice, Polyplax asiatica and Polyplax spinulosa. Minichromosomes are separated by a comma (,). Genes on the same minichromosome 
are linked by a hyphen (-). Genes underlined have an opposite orientation of transcription relative to those not underlined. Minichromosomes 
and gene arrangements that are ancestral to the Polyplax rat lice, the human lice and the pig lice are indicated with solid-line boxes. tRNA genes 
that were translocated in the Polyplax lice are indicated with dot-line boxes. The two minichromosomes of Po. asiatica where the tRNA-leucine 
genes swapped identities are also indicated with dot-line boxes. The phylogeny and the date estimate for the common ancestors are after Light et al. [3]. 



the swap of identities between the two tRNA-leucine genes 
in Po. spinulosa after the split. 

RrnS-trnC minichromosome is present in the rat 
louse, Po. asiatica, and the two Haematopinus species of 
the pig lice, indicating this minichromosome is ancestral 
to the rat lice, the pig lice and the human lice (Figure 4). 
In the rat louse, Po. spinulosa, however, trnS r trnS 2 has 
been inserted upstream rrnS gene (Figure 2B; Figure 4). 
The gene arrangement, nadl-trnT-trnG-nad3-trnW (genes 
in bold have opposite orientation of transcription to others), 
is present in Po. asiatica and the Haematopinus pig lice, 
indicating this arrangement to be ancestral to the rat lice, 
the pig lice and the human lice (Figure 4). In Po. spinulosa, 
however, trnT has been inserted upstream trnD gene, 
leaving a 54-bp pseudo-£raT between nadl and trnG 
(Figure 2B; Figure 5A). The gene arrangement, cox3-trnA, 
which is present in Po. asiatica, the Haematopinus pig 
lice and the human lice, is apparently ancestral to 
these lice. In Po. spinulosa, however, trnA has been inserted 
downstream nad6 gene (Figure 2B; Figure 4). 

Translocation of tRNA genes between minichromo- 
somes also occurred in Po. asiatica. The gene arrangement, 
trnP-cox3, which is present in Po. spinulosa and the 
Haematopinus pig lice, is apparently ancestral to the 



rat lice, the pig lice and the human lice (Figure 4). In 
Po. asiatica, however, trnP has been inserted downstream 
trnN, leaving a 57-bp pseudo-^mP between nad4L and cox3 
(Figure 2 A; Figure 5B). Unusually, there is a 189-bp non- 
coding region in between trnN and trnP (Figure 2A). The 
gene arrangement, nadS-trnF, which is ancestral to insects 
[12], is retained in Po. spinulosa and the Haematopinus pig 
lice. In Po. asiatica, however, trnF has been inserted 
downstream nad4 gene (Figure 2 A; Figure 4). 

The rat louse, Po. asiatica, and the Haematopinus pig lice 
retained an ancestral gene arrangement of insects, coxl-trnL 2 
(Figure 4) [12]. The gene arrangement, trnL r rrnL, which is 
present in Po. asiatica, the Haematopinus pig lice, and the 
human head and body lice, is apparently ancestral to these 
lice. In Po. spinulosa, however, these two trnL genes swapped 
their positions: trnL 1 lies downstream coxl in one minichro- 
mosome whereas trnL 2 lies upstream rrnL in another mini- 
chromosome (Figure 2B). Given the high sequence similarity 
between trnh 1 and trnL 2 (Figure 6A; also see below), it is 
more likely that these two tRNA genes swapped their 
identities by recombination or point mutations at the third 
anti-codon positions than by translocations of these two 
genes between minichromosomes, as observed previously 
in the human head louse and body louse [11]. 



(A) 



(68 bp, 49% similarity) 



PstrnT AGTCCTCACAGTCTTTGATGACTGAGAGTTTGTACCTCTCAAGAGGTTTAGCCTGAGG- -ACTA 62 

Ps pseudo-trnT CGCCTATTAA GC-TCCCTAGCGCCTTTGAGGCTGAGACAAA TGTCATTAGGCA 54 



(B) 



(64 bp, 38% similarity) 



PatrnP CAGGGGGTAATTGATATACAATTTTAACTTTGGGAGTTAAAA GACTCCACTTGGAGTCCCCTGA 64 

Pa pseudo-trnP - - AA GAGAGGGCGTCAGGGATTTCGGAG- - - GAGGGTTTATAGCCCTCTTATAGCCTAGGCT- - 57 

Figure 5 Alignment of the nucleotide sequences between trnT and pseudo-frnT of Polyplax spinulosa (A) and between trnP and 
pseudo-frnP of Polyplax asiatica (B). 



Dong et al. BMC Genomics 2014, 15:44 
http://www.biomedcentral.eom/1 471 -21 64/1 5/44 



Page 8 of 12 



(A)d 




lb- X 'ii 



s i s 



8 II S 



.v'pj 5 -^'p ^/p o=f?^ 



U S J ""J : ~ s 



/If r- s s"s* ft? »r- ! s"s f« s f*f r-° 

" ¥ " ¥ " If " ¥ " ¥ " r 

"J »cJ L S 8 S 8J 8J 



p f£ p; pj pj f v .SJ ,-w' p; 



S 8 g S J us 



rx" r =:x' f (X" n ! {s; r : XC" P- 5 X F 

r r i i s i 

u J s cJ U U u u 



If B if If w 



8 § S 8 Si 



S 8 X 



i s ' s s s 



J 8 S 



1 ir 

US U A 



(B) 



C 5 



if 



Figure 6 (See legend on next page.) 
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(See figure on previous page.) 

Figure 6 Inferred secondary structures of the mitochondrial tRNAs of Polyplax asiatica (Pa) and Polyplax spinulosa (Ps) (A) and 
mitochondrial tRNAs for alanine and cysteine of the wallaby louse, Heterodoxus macropus (B). Identical sequences shared between the 
two tRNAs for leucine are in bold, so are the identical sequences shared between tRNA-alanine and tRNA-cystine (see also Table 2). The 17-bp 
identical sequences shared between the two tRNAs of Heterodoxus macropus are in bold (see also Table 2). 



Recombination between mt genes in parasitic lice 

In the human lice and the pig lice, ten pairs of mt genes 
share stretches of identical sequences much longer than 
expected by chance, providing striking evidence for 
DNA recombination between mt genes and between 
minichromosomes [10-12]. As in the human lice and the 
pig lice, trnLj and trnL 2 in both rat lice, Po. asiatica and 
Po. spinulosa, share 28-bp and 25-bp identical sequences 
respectively (plus a 11 -bp identical sequence immediately 
downstream with 1 bp in between in both species), which 
are 4-5 times longer than in the animals that have the 
typical mt genome organization (Table 2; Figure 6A). 
Recombination between trnLj and trnh 2 appears to be 
common in the blood-sucking lice and occurs likely 
more frequent than other mt genes. 

Two other pairs of mt genes in the rat louse, Po. 
spinulosa, also share identical sequences, 36 bp between 
cob and nadS, and 32 bp between trnA and trnC, which 
are 3-5 times longer than in other animals (Table 2). 
Recombination between cob and nadS, and between trnA 
and trnC in Po. spinulosa, appears to be recent evolutionary 
events and less common as these two pairs of genes do not 
share longer-than-expected sequences in the other rat 



louse, Po. asiatica, nor in the human lice and the pig lice 
(Table 2). Intriguingly, trnA and trnC of the wallaby louse, 
Heterodoxus macropus, share 17 bp identical sequence 
(Figure 6B; Table 2), which is 2-3 times longer than 
expected by chance, indicating that recombination between 
mt genes may also occur in the animals that have the 
typical mt genome organization although it is much 
rarer than in the blood-sucking lice, which have fragmented 
mt genomes. 

It was suggested that recombination between tRNA 
genes could affect the secondary structure of their 
corresponding tRNAs [11]. This is, indeed, the case 
for tRNA- Alanine and tRNA- Cystine in the rat lice. 
The 32-bp identical sequence shared between trnA 
and trnC in P. spinulosa falls on the AA-arm and T-arm 
of tRNA-Alanine and tRNA-Cystine. These two tRNAs 
have the longest T-arms among all of the tRNAs of Po. 
spinulosa, with 8 pairs and 7 pairs respectively at the 
T-stems and 9 nt at the T-loop, in comparison to the 
average 4.2 pairs and 4.6 nt of the 22 tRNAs (Table 3; 
Figure 6A). In contrast, in the other rat louse, Po. asiatica, 
in which trnA and trnC do not share longer-than-expected 
sequences, the T-arms of tRNA-Alanine and tRNA-Cystine 



Table 2 The longest stretches of identical sequence shared by mitochondrial genes in two rat lice, three human lice, 
two pig lice, which have fragmented mitochondrial genomes, and six other species of bilateral animals that have the 
typical mitochondrial genomes 



Pair of gene 



The longest stretches of identical sequences shared 







Rat lice 




Human lice 




Pig lice 


Animals with typical mt genome organization 


Pa 


Ps 


Pc 


Ph 


Pp 


Has 


Haa 


Bm 


Cb 


Hm 


Dy 


Ce 


Hos 


trnLj 


trnL 2 


28 


25, 11 


33, 32 


33, 32 


35, 32 


16, 10,9 


16, 10, 9 


7 


6 


7 


10 


6 


6 


cob 


nad5 


12 


36 


12 


12 


13 


13 


13 


12 


16 


14 


13 


13 


12 


trnA 


trnC 


6 


32 


6 


6 


6 


6 


6 


NA 


7 


17 


7 


8 


6 


nad4 


nad5 


12 


18 


127, 30 


127, 30 


NA 


12 


12 


13 


15 


15 


16 


14 


11 


nad5 


rrnL 


13 


13 


99 


99 


10 


11 


10 


12 


14 


13 


15 


16 


10 


trnG 


trnR 


5 


5 


28, 14 


28, 14 


32, 26 


5 


5 


5 


6 


7 


6 


8 


6 


coxl 


nad4L 


9 


11 


10 


10 


29 


11 


11 


13 


11 


14 


13 


12 


10 


nad2 


rrnL 


10 


11 


26 


26 


10 


10 


11 


13 


11 


14 


13 


12 


10 


trnP 


trnT 


6 


7 


7 


7 


NA 


26 


26 


6 


8 


8 


9 


10 


7 


atp8 


trnG 


7 


9 


26 


26 


9 


6 


6 


10 


11 


11 


12 


NA 


6 


otp8 


nad2 


12 


10 


10 


10 


8 


25 


25 


10 


14 


12 


14 


NA 


11 


trnl 


trnT 


10 


8 


6 


6 


16 


6 


6 


6 


5 


7 


7 


9 


6 



Note: Abbreviations of species names are: Pa, Polyplax asiatica (louse of the greater bandicoot rat); Ps, Polyplax spinulosa (louse of the Asian house rat); Pc, Pediculus 
capititis (human head louse); Ph, (human body louse); Pp, Pthirus pubis {human pubic louse); Has, Haematopinus suis (domestic pig louse); Haa, Haematopinus apri 
(wild pig louse); Bm, Bothriometopus macrocnemis (screamer louse); Cb, Campanulotes bidentatus (pigeon louse); Hm, Heterodoxus macropus (wallaby louse); Dy, Drosophila 
yakuba (fruitfly); Ce, Caenorhabditis elegans (roundworm); Hos, Homo sapiens (human); NA, not available. Stretches of shared identical sequences longer than expected by 
chance are in bold. 
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Table 3 Length of T-stem, T-loop and T-arm of the mitochondrial tRNAs of the rat lice, Polyplax asiatica [Pa) and 



Polyplax spinulosa (Ps) 


tRNA 


Anti-codon 


Pa T-stem (pairs) 


Ps T-stem (pairs) 


Pa T-loop (nt) 


Ps T-loop (nt) 


Pa T-arm (nt) 


Ps T-arm (nt) 


Alanine (A) 


UGC 


5 


8 


7 


9 


17 


25 


Cystine (C) 


GCA 


4 


7 


5 


9 


13 


23 


Aspartic acid (D) 


GUC 


4 


4 


3 


4 


11 


12 


Glutamic acid (E) 


UUC 


5 


3 


11 


4 


21 


10 


Phenylalanine (F) 


GAA 


5 


4 


6 


4 


16 


12 


Glycine (G) 


UCC 


5 


3 


5 


6 


15 


12 


Histidine (H) 


GUG 


4 


5 


4 


3 


12 


13 


Isoleucine (1) 


GAU 


5 


3 


6 


4 


16 


10 


Lysine (K) 


UUU 


5 


3 


3 


3 


13 


9 


Leucine (1^) 


UAG 


5 


4 


4 


3 


14 


11 


Leucine (L 2 ) 


UAA 


5 


4 


6 


4 


16 


12 


Methionine (M) 


CAU 


3 


4 


10 


3 


16 


11 


Asparagine (N) 


GUU 


5 


3 


4 


3 


14 


9 


Proline (P) 


UGG 


4 


4 


4 


3 


12 


11 


Glutamine (Q) 


UUG 


5 


5 


6 


6 


16 


16 


Arginine (R) 


UCG 


5 


3 


6 


4 


16 


10 


Serine (SO 


UCU 


5 


5 


8 


10 


18 


20 


Serine (S 2 ) 


UGA 


5 


5 


6 


4 


16 


14 


Threonine (T) 


UGU 


3 


4 


5 


3 


11 


11 


Valine (V) 


UAC 


5 


3 


9 


4 


19 


10 


Tryptophan (W) 


UCA 


6 


5 


8 


4 


20 


14 


Tyrosine (Y) 


GUA 


3 


3 


4 


4 


10 


10 


Total 




101 


92 


130 


101 


332 


285 


Mean 




4.6 


4.2 


5.9 


4.6 


15.1 


13.0 



have 5 pairs and 4 pairs respectively at the T-stems, 
and 7 nt and 5 nt at the T-loop, which are very close to 
the average 4.6 pairs and 5.9 nt of the 22 tRNAs (Table 3). 

A possible link between the extent of mt genome 
fragmentation and the length of life cycle 

The rat lice in the genus Polyplax share their MRCA 
with the human lice -47 Mya, and share their MRCA 
with the pig lice -67 Mya (Figure 4) [3]. The mt genomes 
of the two Polyplax species we sequenced in the present 
study have 3.4 genes on average on each minichromosome 



and thus, are less fragmented than those of the human lice 
(2.1 and 2.4 genes per minichromosome) but are more 
fragmented than those of the pig lice (4.1 genes per 
minichromosome). Clearly, the extent of mitochondrial 
genome fragmentation varies among different lineages of 
the blood-sucking lice. 

The pig lice, Haematopinus suis and Haematopinus 
apri, have the least fragmented mt genomes among the 
blood-sucking lice whose mt genomes have been sequenced 
completely or near completely sequenced to date [12]. 
Intriguingly, these two Haematopinus species have the 



Table 4 Body length, life-cycle length and number of mitochondrial genes per minichromosome of five species of 
blood-sucking lice 

Species of lice Pig louse Rat louse Human pubic louse Human body louse Human head louse 

Haematopinus suis Polyplax spinulosa Pthirus pubis Pediculus humanus Pediculus capitis 

Adult body length (mm) 4-6 0.6-1.5 1.3-2.0 2.3-3.6 2.1-3.3 

Length of life cycle (days) 29-48 25-28 1 6-25 1 4-20 13-17 

Mitochondrial genes per 4.1 3.4 2.4 2.1 2.1 

minichromosome 

Note: The information of body length and life cycle was from Dynum et al. [31], Florence [32], Wall and Shearer [30] and Baker et al. [33] for the pig louse; from 
Baker [22] and Goodwin et al. [34] for the rat louse; from Buxton [35] and Service [36] for the human pubic louse; and from Buxton [35], Buxton [37] and Service 
[36] for the human body louse and head louse. 
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largest body size (up to 6 mm long) among the blood- 
sucking lice known [30]. The body size of the blood- 
sucking lice, however, does not appear to be linked to the 
extent of mt genome fragmentation as the human 
lice, Pediculus humanus and Pediculus capitis, have 
larger body size (up to 3.6 mm long) than the Polyplax rat 
lice (up to 1.5 mm long) but have more fragmented mt 
genomes than the Polyplax rat lice (Table 4). 

Instead, the length of life cycle of the blood-sucking lice 
appears to be linked to the extent of mt genome fragmenta- 
tion. The pig louse, H. suis, has the longest life cycle 
(29-48 days depending on weather; no data available 
to H. apri) and the least fragmented mt genome among the 
blood-sucking lice whose mt genomes have been sequenced 
completely or near completely to date (Table 4). In contrast, 
the human head louse, Pe. capitis, and the body louse, Pe. 
humanus, have the shortest life cycle (13-20 days) and the 
most fragmented mt genomes. The rat louse, Po. spinulosa, 
has both the length of life cycle (25-28 days; no data avail- 
able to Po. asiatica) and the extent of mt genome fragmen- 
tation in between the pig lice and the human head and 
body lice. Furthermore, the human pubic louse, Pt pubis, 
has both the length of life cycle (16-25 days) and the extent 
of mt genome fragmentation in between the rat lice and the 
human head and body lice. The inverse link between the 
length of life cycle and the extent of mt genome fragmenta- 
tion, indicated in the present study by limited number of 
species, can be tested further with data from more species 
of blood-sucking lice, and other eukaryotes that have 
fragmented mt genomes such as the ichthyosporean pro- 
tists [16], diplonemid protists [17] and box jellyfish [18]. 
Currently, there is no sufficient data and information, to 
our knowledge, to link the length of life cycle (or other life- 
history factors) with the extent of mt genome fragmentation 
in these eukaryotes. Theoretically, if mitochondrial genome 
fragmentation was a continuous process and each generation 
contributed approximately equally towards fragmentation, 
then it would be expected that the extent of fragmentation 
be inversely linked to the length of life cycle. 

Conclusions 

We sequenced the mt genomes of Po. asiatica and Po. 
spinulosa, collected from the greater bandicoot rat, B. 
indica, and the Asian house rat, R. tanezumi. We identified 
all of the 37 mt genes common to animals in these rat lice. 
The mt genes of Po. asiatica and Po. spinulosa are on 11 
minichromosomes; each minichromosome is 2-4 kb long 
and has 2-7 genes. The two rat lice share the same pattern 
for the distribution of the protein-coding genes and 
rRNA genes over the minichromosomes, but differ in 
the pattern for the distribution of 8 of the 22 tRNA 
genes. The mt genomes of the two Polyplax rat lice 
have 3.4 genes, on average, on each minichromosome 
and thus, are less fragmented than those of the human lice 



(2.1 and 2.4 genes per minichromosome) but are more frag- 
mented than those of the pig lice (4.1 genes per minichro- 
mosome). Our results revealed distinct patterns of mt 
genome fragmentation between the two Polyplax species of 
rat lice. We also showed that the extent of mt genome frag- 
mentation appears to have an inverse link with the length 
of life cycle of the rat lice, the human lice and the pig lice. 
Whether, or not, the extent of mt genome fragmentation is 
indeed linked to the length of life cycle, however, needs to 
be tested further with data from more species and broader 
phylogenetic ranges of blood-sucking lice, and other 
eukaryotes that have fragmented mt genomes. 

Availability of supporting data 

The nucleotide sequences of the mt genomes of the rat lice 
supporting the results of this article have been deposited 
in GenBank [accession numbers KF647751-KF647772; 
http://www.ncbi.nlm.nih.gov/] . 
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mitochondrial genomes of the rat lice, Polyplax asiatica [Pa) and 
Polyplax spinulosa (Ps). 

Additional file 2: PCR primers used to verify the mitochondrial 
minichromosomes of the rat lice, Polyplax asiatica and Polyplax 
spinulosa. 

Additional file 3: Alignment of nucleotide sequences of parts of the 
non-coding regions upstream (A) and downstream (B) of the coding 
regions of the 1 1 mitochondrial minichromosomes of Polyplax 
asiatica, the louse of the greater bandicoot rat, Bandicota indica. 

57F and 57R are the primers used to amplify the coding regions of all 
mitochondrial minichromosomes of Polyplax asiatica. 

Additional file 4: Alignment of nucleotide sequences of parts of the 
non-coding regions upstream (A) and downstream (B) of the coding 
regions of the 1 1 mitochondrial minichromosomes of Polyplax 
spinulosa, the louse of the Asian house rat, Rattus tanezumi. 301 F 

and 301 R are the primers used to amplify the coding regions of all 
mitochondrial minichromosomes of Polyplax spinulosa. 
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